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Abstract 

We calculate the electromagnetic (EM) form factors of the pseudoscalar mesons 
in the light-front framework. Specifically, these form factors are extracted from the 
relevant matrix elements directly, instead of choosing the Breit frame. The results show 
that the charge radius of the meson is related to both the first and second longitudinal 
momentum square derivatives of the momentum distribution function. The static 
properties of the EM form factors and the heavy quark symmetry of the charge radii 
are checked analytically in the heavy quark limit. In addition, we use the Gaussia-type 
wavefunction to obtain the numerical results. 


1 Introduction 


The understanding of the electromagnetic (EM) properties of hadrons is an important topic, 
and the EM form factors which are calculated using non-perturbative methods are the useful 
tool for this purpose. There have been numerous experimental [1-7] and theoretical studies 
[8-13] of the EM form factors of the light pseudoscalar mesons (n and K). However, dne to 
difficulties in the experiments, the EM form factors of light vector mesons (p and K*) have 
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fewer investigations than their pseudoscalar counterparts 


even though they could 


provide much information about the bound-state dynamics. As for the EM form factors of 
heavy mesons (which containing one heavy quark), there are much fewer studies than the 
light ones. In the heavy hadron investigation, however, the heavy quark symmetry (HQS) 
II is a fundamental and model-independent property. In this work, we will study the EM 
form factors of the light and heavy pseudoscalar mesons in the light-front framework. We 
will also check whether HQS is satisified or not among these EM properties of the heavy 
mesons. 

The light front quark model (LFQM) is the only relativistic quark model in which a 
consistent and fully relativistic treatment of quark spins and the center-of-mass motion 
can be carried out. Thus it has been applied in the past to calculate various form factors 
[16-22], This model has many advantages. For example, the light-front wavefunction is 
manifestly boost invariant as it is expressed in terms of the momentum fraction variables 
“+” component) in analog to the parton distributions in the infinite momentum frame. 
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Moreover, hadron spin can also be relativistically constructed by using the so-called Melosh 
rotation fl24|. The kinematic subgroup of the light-front formalism has the maximum number 


of interaction-free generators including the boost operator which describes the center-of-mass 
motion of the bound state (for a review of the light-front dynamics and light-front QCD, see 


The paper is organized as follows. In Sec. 2, the basic theoretical formalism is given 
and the decay constant and the EM form factors are derived for pseudoscalar mesons. In 
Sec. 3, we take the heavy quark limit to check whether HQS is satisfied or not. In Sec. 4, 
the numerical result are obtained by choosing the Gaussian-type wavefunction. Finally, a 
conclusion is given in Sec. 5. 


2 Framework 

A meson bound state consisting of a quark qi and an antiquark q 2 with a total momentum 
P and spin S can be written as 

| M(P,S,S Z )) = J{d 3 p 1 }{d 3 p 2 }2(2n) 3 S 3 (P-p 1 -p 2 ) 

x y SSz (Pi,P 2 , Ai, A 2 ) ki(pi,Ai)g 2 (P 2 ,A 2 )), (1) 

Ai,A2 

where p\ and p 2 are the on-mass-shell light-front momenta, 

P = {P + ,P±) , P± = ( P l ,P 2 ) , P~ = — + + P± , ( 2 ) 

p+ 


2 












and 


3 _ dp + d 2 pj 


{d P} = 2 ( 2^)3 ' 

k(Pi,Ai)g(p 2 ,A 2 )) = &A! (Pi) d A 2 (Pz) 10), 

{h'(p'),b[(p)} = {d y (p’),d[(p)} = 2(2t r) 3 5 3 {p' -p) 5y\. 

In terms of the light-front relative momentum variables (x, k±) defined by 

pf = ( 1 - x)P + , pj = xP + , 

Pi± = (1 - x)P± + k±, p 2 ± = xP ± - k±, 

the momentum-space wavefunction d' 5,s ' z can be expressed as 

V SS ’(pi,P 2 , Ai, A 2 ) = R\^\ 2 (x, k±) <j)(x, k ± ), 


( 3 ) 


(4) 


(5) 


where 0 (x, k±) describes the momentum distribution of the constituents in the bound state, 
and R\^x 2 constructs a state of definite spin (5, S z ) out of light-front helicity (Ai, A 2 ) eigen¬ 
states. Explicitly, 

R\i\ 2 (x,k±) = J2( X ^ n M( 1 - x , k x, m i)\si)(X2\'RM( x ,- k x,m 2 )\s2)(-s 1 --S2\S,S z ), ( 6 ) 

Sl,S2 ^ ^ 


where |sj) are the usual Pauli spinors, and TZm is the Melosh transformation operator 

rrii + xM 0 + id ■ k± x n 


Rm(x, k ± ,rrii) = - , 

y {mi + xM 0 ) 2 + 

with n = ( 0 , 0 , 1 ), a unit vector in the ^-direction, and 

M 2 = rnj + kj mf+fcj . 

0 (1 — x) x 

In practice, it is more convenient to use the covariant form for R\f\ 2 


(7) 


( 8 ) 


R S \i\ 2 ( x , k ±) = 


pi pi 


x/2 Mq 


u{pi, Ai)rn(p 2 , A 2 ), 


(9) 


where 

M 0 = \/.U ( ? - (m.\ - rth,) 2 . 
r = 75 (pseudoscalar, 5 = 0). 

We normalize the meson state as 

<M(P', 5', 5')|M(P, 5, S z )) = 2(2 tt) 3 P + 8 3 (P' - P)5 s , s 8 s , Sz 


( 10 ) 
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so that the normalization condition of the momentum distribution function can be obtained 

J{dx} \<j)(x, k ± )\ 2 = 1 , ( 11 ) 

where 

, 1 dxd 2 k± 

{dx} = 

In principle, the momentum distribution amplitude <p(x, k ±) can be obtained by solving the 
light-front QCD bound state equation [|5|]. However, before such first-principles solutions are 
available, we would have to be contented with phenomenological amplitudes. One example 
that has been often used in the literature for heavy mesons is the Gaussian-type wavefunction, 


. . r dk- ( k 2 

G= A/V—expU— 


( 12 ) 


where J\f = 4 ( 7 r/cn 2 ) 3//4 and k~ is of the internal momentum k = (k±, k z ), dehned through 

ei - k z 


l-x = 

ei + e 2 

with ei = \/m 2 + k 2 . We then have 

Mq = e\ + e 2 , k z = 

dk z 


x = 


and 


2 

ere 2 


e 2 + k z 

ei + e 2 ’ 

(13) 

+ k\ 

2xM 0 ’ 

(14) 

5 

(15) 


which is the Jacobian of transformation from (x, k ±) to k. 

2.1 Decay Constants 

The decay constant of a pseudoscalar meson P(gig 2 ) is dehned by 

( 0 | 2 l„|P(p)) = ifppp, 


(16) 


where A M is the axial-vector current. It can be evaluated using the light-front wavefunction 
given by (0) 

(0|g 2 7 + 7 5 gi|P) = J{d 3 p 1 }{d 3 p 2 }2(2TT) 3 S 3 (p-p 1 - p 2 )0p(a:, k ± )R° x ° lX2 (x, k±) 

x ( 0 |g 2 7 + 7 5 gi|gig 2 ). (17) 


Since M 0 ^x( 1 — x) = \JA 2 + k'j_, it is straightforward to show that 

Gt [ {dx} ^kM At 


fp = 


V 2 J 


(18) 


where A = m\X + m 2 (l — x). Note that the factor \f?> in (|T£|) arises from the color factor 
implicitly in the meson wavefunction. 
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2.2 Electromagnetic Form Factors 

The EM form factor of a pseudoscalar meson P, Fp(Q 2 ), is determined by the scattering of 
one virtual photon and one meson. It describes the deviation from the point-like structure of 
the mesonch, and is a function of the square of the photon momentum Q. Here we consider 
the momentum of the virtual photon in space-like region, so it is always possible to orient 
the axes in such a manner that Q + = 0. Thus the EM form factor is determined by the 
matrix element 


(P(P')|J + |P(P)> = e F P (Q 2 )(P + P') + , 


(19) 


where = ge^ey^g is the vector current, e q is the charge of quark q in e unit, and Q 2 = 
— (P' — P) 2 > 0. With LFQM, Fp can be extracted by Eq. 


F P (Q 2 ) = e 


91 


e 92 


(p P {x,k ± ) <j) P ,(x,k' ± ) 

sjA 2 + k\ \JA 2 + kf 
(p P (x,k±) (f)p'(x,kl) 


A 2 + k _l • k '_l 
A 2 + k _ L • k'[ 


( 20 ) 


sjA 2 + k 2 ± sj A 2 + k” 2 

where k' L = k±+xQ ±, k’[_ = k± — (1 — x)Q±. From Eqs. (i. ©■ and (|^), it is understandable 
that the term \JA 2 + k]_ comes from the Melosh transformation. After fixing the parameters 
which appear in the wavefunction, Eq. (j20|) can be used to fit the experimental data. But 
this is not the whole story. We consider the term 0p = (f>p(x, k±)/^JA 2 + k 2 L and take the 
Tayor expansion around kp 


(,i>p’{k'l) = <j> P >(kj_) + 


deppr 


dki 


(k't - kl) + 


d 2 <pp / 


Qx=o 


2(dFlf 


0 k’l - kl) 2 + 


Q ±=o 


Then, by using the idenity 

Id 2 k ± Ik, ■ A i )!k - D ) = ) jd 2 k ± kl A ± ■ B i. 


( 21 ) 


( 22 ) 


we can rewrite 
->2 


as 


Fp(Q ) — ( e qi +e q2 ) 


91 1 1 ' 92 ) 

Q 2 [{dx} (j) 2 P (x, k_ l) [x 2 e qi + (1 - x) 2 e q2 ] f©p ^ n + + ®pk 2 


A 2 + k 2 


where 


0(Q 4 ), 


©M — 


1 f d(p M ^ _ 1 f d 2 (j) M 


<t>M \ dkl J 


4>m V (dk 2 ± y 


(23) 


(24) 


From Eq. (|23|), the static property Pp(0) = ep is quite easily checked. The mean square 
radius of the meson P is determined from the slope of Fp at Q 2 = 0: 

(rV ,_ 6 ^« 2 ) 


dQ 2 


(25) 


Q 2 =0 
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It should be realized that the size and the density of a hadron depend on the probe. For an 
electromagnetic probe, it is the electric charge radius that is obtained. From the experimental 
view, (r 2 )p cannot be measured directly and is obtained by fitting the data on Fp to a pole 
or dipole form. Here we easily obtained the equation of (r 2 )p 


(r 2 ) P = (r 2 ) qi + (r% 


92 


= -6 


+ e q A -6 


X 2 (j)p 




— x) 2 4>p {A~ + + (A 2 + k 2 ± )k 2 ± {-^2 


.(26) 


From Eq. (p 6 |), it is worthwhile to mention that, hrst, the mean square radius of a me¬ 
son is the sum of the contributions of the valence quarks. Second, (r 2 ) is related to the 
hrst and second longitudinal momentum square derivatives of 0 which contain the Melosh 
transformation effect. 


3 Heavy Quark Limit 


In this section, we will check the HQS among the charge radii by taking the heavy quark limit. 
To proceed, we hrst investigate the heavy-quark-limit behavior of the wavefunction. Since 
the x in the normalization condition ( 11 ]) is the longitudinal momentum fraction carried by 
the light antiquark, the meson wavefunction should be sharply peaked near x Aqcd/ mg' 
It is thus clear that only terms of the form “mgo;” survive in the wavefunction as mg —> oo; 
that is, rriQX is independent of itiq in the itiq —> oo limit. In the itiq —> oo limit, we must 
rewrite Eq. (|T0|) in the mQ-independent form 



d 2 k _l 
2(2vr) 3 


|$( X ,^)| 2 = 1 , 


where X = rriQX and |26 


Hx, k ± ) 




(27) 


(28) 


The scaling behavior of Eq. (^) is the constraint of the light-front wavefunction when we 
consider the infinite quark mass limit. For the Gaussian-type wavefunction (0. it satisfies 
an asymptotic form 


®(X, k A 




3/4 




exp 


2w 2 


exp 



(29) 


Thus we can use this wavefunction when the heavy quark limit is considered. 

In the rriM, tuq —> oo limit it is appropriate to describe the meson state with the meson 
velocity v 


H 


I «W> = ">m 1/2 | M(P)), 


(30) 
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where v = P/rriM- For the decay constant, the definition (|I 6 D becomes 

(0|g7M75Q|^O)) = i fp v n: (31) 

and in the mg —> oo limit it is 

A 


- a/3 /• dX d 2 k± 

fp = 4 —= / - 7 -<F(A, k A 

y/2J 2 27r) 3 V ’ 


A 2 + k 2 


(32) 


where M = X + m g - 2 . Comparing Eq. (p?2| ) with Eq. (|T%|), we obtain the HQS scaling law for 
the decay constant: 


fp — 's/ppm fp- (33) 

For the mean square radius Eqs. (|26|) , when the heavy quark limit is considered, we obtain 

(r 2 )p=(r 2 ) Q + (r 2 ) g - 2 , (34) 

where 


( r )q = e Q\ZP2 


-6 r dXd 2 k± 
m Q J 2 ( 2 tt ) 3 


X 2 $ (^4 2 + 2 k 2 ± )J^ + (^4 2 + kl)k] 


d 

dff 




$ 


0, 


(35) 


[r 2 )q 2 = eg 2 { -6 


dXd 2 k A 

2(2vr) 3 


■ 4 > 


{A 2 P2kl)4r + ^ 2j rkl)k 2 A 
dk~ t 


dk 2 . 


$ 


(36) 


and <F = <F/ \JA 2 + k 2 ± . Eq. (|35D means that the mean square radius ( r 2 )p is blind to 
the flavor of Q. This is the so-called flavor symmetry. We find that the light degrees of 
freedom are blind to the flavor of the heavy quark. In addition, Ref. |27] finds the mean 


square radius also satisfied the spin symmetry. These are the so-called HQS. Up to now, we 
have not used the wavefunction yet, this also satisfies the well-known property that HQS 
is model-independent. Reviewing the processes, we can realize that, in this approach, the 
static properties of the EM form factors and the heavy quark symmetry of the mean square 
radii can be checked much more easily than in the Breit frame. This is the major reason 
why we calculate the Q 2 dependence of those form factors order by order. 

We must emphasized here that, in the mQ —> oo limit, the vanishing of the heavy quark 
sector in the form factor is true only for the Q 2 —> 0 region. In the time-like region, near the 
threshold for the meson pair production the heavy quark sector is dominant and described 
by the Isgur-Wise function |28| . 


4 Numerical Results 


In this section, we will use the Gaussian-type wavefunction ( [12]) to calculate the EM form 
factors and the mean square radius. The parameters appearing in the wavefunction, the 
quark mass m q and the scale parameter u, are constrained by the decay constants. 
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The decay constants of the pseudoscalar mesons n and K come from experiments |29 

f n = 130.7 MeV, f K = 159.8 MeV, 
the others are obtained by lattice and constituent quark model: 


( 37 ) 


f D = 192 MeV[30], f D , = 210 MeV[g, f„ = 157 MeVgO 
fa. = 171MeVBg, } B = 360 MeVED. 


(38) 


Combining with the quark masses 

m Ut d = 0.24 GeV, m s — m Ut d = 0.18 GeV, m c = 1.6 GeV, nib = 4.8 GeV, (39) 
we fit the scale parameters 


= 0.333 GeV, u K = 0.379 GeV, u D = 0.443 GeV, u Ds = 0.450 GeV, 
co B = 0.477 GeV, uj Bs = 0.485 GeV, u Bc = 0.813 GeV. (40) 


There are differences between these parameters and the ones in |3T[ because the wavefunc- 
tions in the two cases are not the same. However, they have a common tendency such that 
u Mi < ojM j if Mi < Mj. This corresponds to the ordering law for the size of heavy-light 


bound states. 

The Q 2 -dependences of F n and F K can be obtained by Eq. (|2D|), and we compare the 
results with the data in Fig. 1 and 2, respectively. In addition, the mean square radii of the 
pseudoscalar meson can be obtained by Eq. (PU|). We list the results of the {r 2 ) n + ,k° an d 

the experimental data in Table 1. (the unit is fin 2 ) 


{r 2 ) 

7T + 

K+ 

K° 

this work 

0.443 

0.349 

-0.0676 

[11] 

0.314 

0.240 

-0.020 

[13] 

0.452 

0.38 

0.057 

experiment 

0.439 ± 0.008 [1] 

0.34 ± 0.05 [5] 

-0.054 ± 0.026 [7] 


Table 1. The mean square radii of the 7r + , K + , and K° mesons. 


The negative signs in Table 1 are interesting, and may be interpreted as the preponderance of 
negative electric charge in the tail of the distribution. We find these values are all consistent 
with the data. Comparing with Ref. [|J, they also used the light-front approach. There were 
various parameter combinations to fit the data of F n for both small and large momentum 
transfers. 

According to the vector meson dominance (VMD) model [|T!|, there is a physical explana¬ 
tion: the pion form factor is determined by a p-meson pole. Generally speaking, this simple 
picture fits the data well. A detailed study |13[ obtained a better fit when one considers the 
p-u mixing and three vector meson (p, u, and 0) poles to the pion and kaon form factors, 
respectively. 

On the other hand, the mean square radii of the heavy pseudoscalar meson have not 
been measured yet. For comparsion, here we define and calculate them as ( t 2 )fm f° r the 
























finite quark masses and as ( r2 )lM for the infinite quark masses. In the case of the infinite 
quark masses, the decay constant fp cannot be measured in the true world, so we obtain it 
approximately by using the values fp = 157 MeV and rnp = 5.28 GeV in Eq. (|33|) . The 
results are listed in Table 2. 



D + 


Dt 

B + 

R u 


Bt 

( T 2 )fm 

0.184 

-0.304 

0.124 

0.378 

-0.187 

-0.119 

0.0433 

(r 2 ) IM 

0.248 

-0.496 

0.181 

0.496 

-0.248 

-0.181 



Table 2. The mean square radii of the heavy pseudoscalar mesons for the finite 
quark masses (' t2 )fm and for the infinite quark masses (r 2 )iM- 


From Table 2, we cannot obviously find the situation that, comparing with the values in 
D q system, the ones in the B q system are closer to those in the infinite-quark-mass system. 
The reason is that the (r 2 ) is sensitive to the fp, but the uncertainty of the decay constant 
In fact, if we use the most recent value fp> s = 280 MeV 0 


is not small. 


the result 
2 \ 


(r 2 ) D + = 0.083 fm" is quite different from the one in Table 2. For the B c meson, the (r 2 )/M 
are not given here because both b and c quarks are heavy. The HQS must be reconsidered 
in this case. 


5 Conclusion 

We have calculated the EM form factors of the pseudoscalar mesons. The EM form factors are 
extracted from the relevant matrix elements directly, instead of choosing the Breit frame. We 
found that the charge radius is related to both the first and second longitudinal momentum 
square derivatives of the momentum distribution function. We also found that the static 
properties of the EM form factors and the heavy flavor symmetry of the mean square radii 
are checked analytically by evaluating the Q 2 dependence of those form factors order by order. 
Therefore, in the heavy quark limit, the charge radii of pseudoscalar have flavor symmetries, 
and these properties are model-independent. In addition, The Q 2 -dependences of the form 
factors F nK and the mean square radius of light and heavy mesons have been calculated 
by using the Gaussian-type wavefunction. The form factors F n and F K in small momentum 
transfer and the values of {r 2 ) n + k +are all consistent with the current experimental data. 
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FIGURE CAPTIONS 


Fig. 1 The charge form factor of the pion in small momentum transfer. Data are taken 
from HU. 

Fig. 1 The charge form factor of the Kaon in small momentum transfer. Data are taken 
from ||. 
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